Crystal field splitting and correlation effect on the electronic structure of v4 2 Ir03 
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The electronic structure of the honeycomb lattice iridates Na2lrC>3 and Li2lr03 has been in- 
vestigated using resonant inelastic x-ray scattering (RIXS). Crystal- field split d-d excitations are 
resolved in the high-resolution RIXS spectra. In particular, the splitting due to non-cubic crystal 
fields, derived from the splitting of j e g=3/2 states, is much smaller than the typical spin-orbit en- 
ergy scale in iridates, validating the applicability of j e s physics in ^blrOs. We also find excitonic 
enhancement of the particle-hole excitation gap around 0.4 eV, indicating that the nearest-neighbor 
Coulomb interaction could be large. These findings suggest that both Na2lr03 and Li2lr03 can be 
described as spin-orbit Mott insulators, similar to the square lattice iridate S^IrCU. 

PACS numbers: 75.10.Jm, 75.25.Dk, 71.70.Ej,78.70.Ck 



The intense interest in iridium oxides, or iridates, arises 
from a number of competing interactions of similar mag- 
nitude [H-Q. While the on-site Coulomb interaction is 
the dominant energy scale in 3d transition metal oxides, 
the spin-orbit coupling (SOC) is largely ignored. On the 
other hand, for 5d elements such as Ir, the SOC becomes 
significant, and in fact plays a dominant role. A good 
example is Sr2ix04, whose electronic states are well de- 
scribed by Jeff =1/2 states arising from the spin-orbit split 
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tig levels 

One of the most intensely scrutinized families of iri- 
dates is the honeycomb lattice family A2I1O3 (A=Na,Li) 
@,0j.l3"13|- Originally thought of as topological insula- 
tor [5j , these materials are now believed to be Mott insu- 
A recent calculation though suggests that 
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uniaxial strain might still drive the system to topologi- 
cal insulating behavior |14j ]. Furthermore, these materi- 



als could be described with the Kitaev-Heisenberg model 
@, E3 j in which bond-dependent Kitaev interaction are 
realized and support various types of topological phases. 
The applicability of such intriguing theoretical possibili- 
ties to real system crucially depends on the j e ff physics 
arising from strong SOC. However, the experimental sit- 
uation seems to be far from clear. In particular, struc- 
tural refinements find a sizable trigonal distortion of the 
IrOg octahedra 15, 16| , which will produce crystal field 
splittings within the t 2g manifold. If the splitting is com- 
parable to the SOC, the j e g=l/2 states will mix with 
jeff =3/2 states and the relevant microscopic model be- 



comes quite different from the ideal j ff physics [13|, [17| , 
preventin g th e Kitaev-Heisenberg model from being re- 
alized M, nil, lia . Recent theoretical studies have even 



suggested that the ground state has a large contribution 
from the j e ff=3/2 state [lij ]. 

Therefore, it is of great importance to elucidate the un- 
derlying electronic structure of Na2lr03 experimentally. 
In particular, the spectroscopic investigation of excita- 
tions between spin-orbit split j e ff states can provide us 
with direct information regarding the size of the crys- 
tal field splitting with respect to the typical SOC energy 
scale in iridates (0.4-0.5 eV). In the case of Si^IrGj, such 
excitations between j e ff=3/2 to j cff =l/2 were observed 
around 0.6-0.8 eV in the resonant inelastic x-ray scatter- 
ing (RIXS) data Q, which is accounted for in the quan- 
tum chemical calculation by Katukuri ct al. [19j . The 
splitting within these "spin-orbit" excitations arises due 
to non-zero tetragonal crystal fields, and is much smaller 
(~ 0.1 eV) than the SOC, justifying the j e ff description 
of Sr 2 Ir0 4 . 

In this Letter, we present a comprehensive picture 
of the low energy electronic structure of Na2lr03 and 
Li2ix03, based on Ir L^-edge RIXS experiments. Our 
high-resolution RIXS measurements allow us to resolve 
the crystal field splitting of the j e ff=3/2 states due to the 
trigonal distortion, which is determined to be about 110 
meV in both compounds. This energy scale agrees very 
well with quantum chemical calculations, and is much 
smaller than the typical value for SOC, validating the 
jeff picture in these compounds. We have also studied 
momentum dependence of the insulating gap; the ob- 
served flat dispersion of the insulating gap is consistent 
with what is expected from a significant Coulomb inter- 
action in both compounds. Taken together, we argue 
that just as Sr 2 Ir04, the honeycomb ^4 2 Ir03 iridates can 
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be described as spin-orbit Mott insulators @, Hi, 17, 2(J 



TABLE I: RIXS and MRCI+SOC excitation energies (C2/m 
structure) for 213 iridates (eV). 
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FIG. 1: (Color online) Top: Wide energy range RIXS spec- 
trum for a single-crystal sample of Na2lrC>3 at Q = (0 6.7) 
obtained with Ei = 11.217 keV. Note the different scale used 
for left and right panel. Bottom: RIXS spectrum for Li2lrC>3 
powder sample at |Q| ~ 8 A~\ obtained with the same Ei. 
All spectra were measured at room temperature. The black 
dashed curves are the result of a fit (see text) and the red 
solid lines represent the background. 

The RIXS experiment was carried out at the Advanced 
Photon Source using the 30ID MERIX and 9ID RIXS 
spectrometer. A spherical (1 m radius) diced Si(844) ana- 
lyzer and Si(844) secondary monochromator were used to 
obtain overall energy resolution (FWHM) of ~ 35 meV. 
To minimize the elastic background, most of the measure- 
ments were carried out in a horizontal scattering geome- 
try near Q = (0 6.7), for which the scattering angle 29 
was close to 90°. We use the C2/m notation for the lat- 
tice jTH , L6 1 . A single crystal of Na2ir03 and a polycrys- 
tallinc samples of Li2lr03wcre grown by the solid state 
synthesis method, previously described in detail [HI EJ ■ 
The Na2lrC>3 crystal was plate-like with a flat shiny sur- 
face, the surface normal was in the (001) direction. 

The RIXS process at the L3-edge of Ir (or any other d 
electron system) is a second order process consisting of 
two dipole transitions (2p — > 5d followed by 5c? — > 2p). 
Therefore, it is especially valuable for detecting excita- 
tions between the d-levels and has been extensively uti- 
lized in the study of 3c? transition metal compounds plT - 

pos- 

n. 



25] . Recent instrumental advances have made it 
sible to measure collective magnetic excitations [1 
In ^2^03, Ir 4+ ions are in the 5c? 5 configuration in a 
slightly distorted octahedral environment of oxygen ions, 
with the edge-sharing IrC>6 octahedra forming a honey- 
comb net. Due to the octahedral crystal field, there ex- 
ists a fairly large splitting (lODq) between the ti g and e g 
states. Since the 5d orbitals are spatially more extended 
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than the 3d orbitals, the lQDq value is expected to be 
much larger. Indeed, in our RIXS investigations of vari- 
ous iridium compounds, well separated tig and e g states 
have been observed, with the 10-Dg value typically about 
3 eV [13. 

In Fig. [1] a representative high-resolution RIXS spec- 
trum of Na2lrC>3 is plotted on a wide energy scale. This 
scan was obtained at room temperature and plotted as 
a function of energy loss {fun = Ei — Et). The incident 
energy, Ei = 11.217 keV, was chosen to maximize the 
resonant enhancement of the spectral features of inter- 
est below 1 eV. A broad and strong feature is observed 
at 2-4 eV and other sharper features are observed be- 
low 1 eV, corresponding to d-d transitions from occupied 
tig states into the empty e g and ti g levels, respectively. 
Also plotted in the figure is the room temperature data 
of polycrystalline Li2lrC>3. Lack of significant momen- 
tum dependence of these d-d excitations (shown later 
in Fig. 3) allows one to directly compare the peak po- 
sitions between the single crystal and powder samples. 
The spectra were fit to 5 peaks (labeled A-E), as shown 
by the black dashed lines. The low energy excitations 
can be fit to three peaks, two Gaussians (B and C) of the 
same width and one Lorentzian (A) on top of a broad 
background (Gaussian). Two Lorentzian functions with 
sloping background were used to fit the higher energy 
excitations (D and E). The resulting peak positions are 
listed in Table Q] 

To clarify the nature of the excitations revealed 
by RIXS, we have carried out multiconfiguration 
self-consistcnt-ficld and multireference configuration- 
interaction (MRCI) calculations 28| on clusters consist- 
ing of one central IrOg octahedron, all adjacent Na or Li 
ions, and the three nearest- neighbor (NN) IrOe octahe- 
dra (see Ref. [H and Supplemental Material for details) . 
Local d-d transitions are computed only for the central 
IrOe octahedron while the NN octahedra are explicitly 
included in the cluster for providing an accurate descrip- 
tion of the nearby charge distribution. Two different lat- 
tice configurations are considered, i.e., the C2/c structure 
[ToL |30| and also the C2/m arrangement proposed more 
recently [HEEIl!- 

Results of spin-orbit MRCI (MRCI+SOC) calculations 
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using the C2/m configuration are listed for Na2lr03 
in the third column of Table I. The MRCI+SOC data fit 
the experiment reasonably well, with peaks B and C cor- 
responding to j c ff =3/2 to j e s=l/2 electronic transitions. 
Above 2.5 eV, the MRCI+SOC results indicate multi- 
ple tig to e g excitations displaying a two-peak structure 
reminiscent of the D and E features in the RIXS spectra. 
However, MRCI+SOC seems to overestimate somewhat 
the relative energies of those latter features. Interest- 
ingly, for the alternative C2/c structure of Na2lr03 [Ioj |. 
the splitting between the two doublets originating from 
the jeff =3/2 quartet in an ideal octahedral environment 
is much larger and the position of the C peak is overes- 
timated by 0.25 eV in the MRCI+SOC treatment. Since 
the deviations from the experimental data are in this 
case larger, the MRCI+SOC results for C2/c symmetry 
are not listed in Table I. The ti g splittings in calculations 
with no SOC are in fact as large as 0.6 eV for the C2/c 
structure of Na2ir03, which gives rise to a highly uneven 
admixture of ti g components in the spin-orbit calcula- 
tions. In contrast, for the C2/m configuration, the ti g 
splittings are about 0.1 eV and the three different ti g 
hole configurations contribute with similar weight to the 
spin-orbit ground-state wave function (see Table II). 

For Li 2 Ir03, the calculations correctly reproduce the 
shift to higher energies of the ti g to e g transitions relative 
to those in Na2lr03. The discrepancy between the ex- 
perimental values and the MRCI+SOC results (e.g., peak 
C) could be caused by the uncertainty in the structural 
model used for this calculation (C2/m from Ref. [29j]). 
Since local structural disorder is not easily captured in 
the regular diffraction data, local structure probes such 
as pair-distribution function (PDF) measurements can 
sometimes be useful for clarifying the structural details. 
We have carried out X-ray PDF studies on Li2lr03 and 
Na2lr03 powder samples. Details of these measurements 
and the comparison of the two structures are reported in 
the Supplemental Material. Except for the overall lattice 
contraction, the Li 2 Ir03 PDF seems to be well described 
by the C2/m symmetry, eliminating the local structural 
disorder as a possible explanation. Most likely cause of 
the structural uncertainty is the oxygen position, since 
x-ray structural probes are not particularly sensitive to 
light elements like oxygen [29|, We note that the 

latest refinements using both powder neutron and single 
crystal x-ray data on Na2lr03 do show important differ- 
ences compared to earlier x-ray powder diffraction data 
and the MRCI+SOC results arc very different for the two 
structures. Better structural refinements using neutron 
diffraction would reduce the oxygen position uncertainty 
in Li2lr03 and could improve the agreement between our 
MRCI+SOC calculation and the experiment. 

One of our main findings is that the splitting of the 
strong RIXS peak located at 0.7-0.8 eV is due to the 
trigonal distortion which is well corroborated with our 
MRCI+SOC calculations. The fact that this splitting 



TABLE II: Percentage contributions of the different Ir 5d 
configurations to the lowest on-site d-d excited states in 
Na2lr03, as obtained from MRCI+SOC calculations. 
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(110 mcV) is much smaller than a SOC of 0.4-0.5 eV 
strongly supports that these excitations are transitions 
from crystal-field-split j e ff=3/2 levels to the j e s=l/2 
state (labeled spin-orbit exciton in Ref. Q). Given that 
the optical gap in this material is about 350 meV [20[ 
and that there is no such excitation in the MRCI+SOC 
calculations which only look at on-site d-d excitations, it 
is reasonable to associate feature A at low energy as aris- 
ing from the excitation of a particle and hole pair across 
the charge gap. Additional periodic density functional 
theory (DFT) calculations shows that a moderate size 
U and SOC can indeed open a (Mott) gap of 300-400 
meV, in accordance with the experimental observation 
(see Supplemental Material). 

The nature of the charge excitation gap can be fur- 
ther revealed by its momentum dependence. In Fig. [21 
we plot the momentum dependence of the low energy 
peaks (A-C) in Na2ir03. In the honeycomb plane, the 
magnetic ordering doubles the unit cell (3l| . and corre- 
spondingly the first Brillouin zone (BZ) becomes smaller. 
Two different BZ schemes are illustrated in the inset of 
Fig. [2] (a) to aid the comparison. We will use the rect- 
angular BZ notation. Note that the two high symmetry 
directions of interest, the q=(H 0) and q=(0 K) in rect- 
angular notation, correspond to the T — M and T — K 
directions in the honeycomb plane, respectively. One can 
see that the overall momentum dependence of the peak 
positions is very small, except for peak A. To investigate 
the behavior of peak A in detail, the low energy portion 
of the spectra was fit to a Lorcntzian peak. Since the 
peak seems to disappear at q=(l 0), we have used the 
spectrum at this q as an empirical background. The fit- 
ting results for peak positions, widths, and intensities are 
shown in Fig. [2fb)-(c) . The width and peak position re- 
mains almost unchanged (sa 10 meV dispersion), but the 
intensity is strongly peaked around the BZ center. This 
can be clearly seen in the pseudocolor plot of the spectra 
shown in Fig. [2jd), in which a strong peak around q=(0 
0) and 0.42 eV is contrasted with the q-independent fea- 
tures B+C. In addition, one can see that the spectral 
weight changes abruptly around 0.4 eV, confirming that 
this is the particle-hole continuum boundary. Based on 
our RIXS results, the electronic excitations in AihOj, 
can be summarized as shown in Fig.[2je). 

It is clear from this observation that the insulating 
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FIG. 2: (Color online) (a) Momentum dependence of the low 
energy RIXS spectra of Na2ir03 obtained at T = 9 K. The 
inset shows a schematic diagram of the (H K 0) reciprocal 
space plane. The Brillouin zones (BZ) corresponding to the 
monoclinic unit cell are blue rectangles. For comparison, we 
also plot the BZ of the honeycomb net in black. The circles 
are the points where RIXS spectra are taken. The low en- 
ergy peaks denoted with red triangles are fit to a Lorentzian, 
and the momentum dependence of (b) the peak position and 
width, and (c) the peak intensity are shown, (d) Same data 
are plotted in false color scale, (e) Schematics of electronic 
excitations in ^IrOa determined from our RIXS measure- 
ments. 



gap is direct (minimum gap at T). The relatively fiat 
dispersion observed in our data is also consistent with 
the DFT calculation which suggests that the correlation 
effect makes the bandwidth smaller, leading to an al- 
most dispersionlcss charge gap. The sharpness in en- 
ergy and momentum of peak A is quite reminiscent of 
the excitonic behavior of the BZ center particle-hole ex- 
citation across the charge-transfer gap in the insulat- 



ing cuprate La2Cu04 [32j. This suggests that an extra 



nearest-neighbor Coulomb interaction V (in addition to 
the on-site interaction U) might be important for mod- 
elling this material. Sizable V could promote the ten- 
dency towards exciton binding and also further narrow 
the bandwidths. The smaller intensity of the charge gap 



feature in Li 2 Ir03 compared to Na2lrC>3 could be due 
to the fact that the Li 2 Ir03 data are powder averaged. 
However, one cannot rule out the possibility of weaker V 
in Li 2 Ir03 as compared to Na 2 Ir03. 

Another interesting aspect of our data is that the dis- 
persion of the gap appears to follow the underlying hon- 
eycomb lattice rather than the crystallographic/magnctic 
unit cell. This is clearly observed by the spectrum ob- 
tained at q=(2 0). While (2 0) is the next BZ center along 
the T—K direction, (1 0) is on the zone boundary; peak A 
disappears at (1 0) but recovers its intensity at the q=(2 
0) position. Additional momentum dependence data, re- 
ported in the Supplemental Material, shows the lack of 
momentum dependence along the L-direction (perpen- 
dicular to the honeycomb plane). Therefore, the elec- 
tronic structure of Na 2 Ir03 seems to be quite well de- 
scribed as that of a 2D honeycomb lattice. 

It is worth comparing the observed low energy RIXS 
spectrum with that of Sr 2 Ir04. In Sr 2 Ir04, a low en- 
ergy magnon was observed below 200 mcV, while highly 
dispersive excitations were observed between 0.4 eV and 
0.8 eV. This latter band of excitations is composed of 
particle-hole excitation across the Mott gap and spin- 
orbit excitations from j e ff=3/2 states to the j e ff=l/2 
states. Because of the smaller single-particle band width 
in j4 2 Ir03 (see DFT calculations in Rcf. Il4j) , the "j e ff 



excitation" in Na 2 IrC>3 is almost dispersionless, unlike 
the highly dispersive counterpart in Sr 2 Ir04. Perhaps 
an even more significant difference is the well separated 
energy scale of the j e g excitation and the particle-hole 
continuum in Na 2 Ir03. These two energy scales are very 
similar in Sr 2 Ir04, but the large separation in Na 2 Ir03 
allows one to investigate these two types of excitations 
separately. 

To summarize, we have carried out a resonant inelastic 
x-ray scattering investigations of electronic excitations in 
Na 2 Ir03 and Li 2 Ir03. Wc observe three well-defined fea- 
tures below 1 eV and a broad two peak feature at 2-5 eV. 
By comparing our observation with quantum chemical 
and density functional theory calculations, we associate 
these features with d-d transitions. Specifically, the high 
energy excitations are from £ 2£( to e g excitations, while 
the low energy excitations around 0.7-0.8 eV are excita- 
tions from j eff =3/2 to j e ff=l/2 states. The splitting of 
the latter feature arising from the trigonal crystal field 
is about 110 meV, much smaller than the spin-orbit cou- 
pling energy scale of Ir compounds, which validates the 
applicability of j e ff physics in j4 2 Ir03. In addition, wc 
observe a lower energy excitation around 0.4 eV, which 
shows very little momentum dependence and is associ- 
ated with the particle-hole excitation across the Mott 
gap; the "excitonic" behavior of this peak suggests the 
nearest-neighbor Coulomb interaction V is sizable. We 
conclude that the electronic structures of both Na 2 Ir03 
and Li 2 Ir03 are similar and these systems can be de- 
scribed as spin-orbit Mott insulators. 



5 



We would like to thank Y. B. Kim and S. Bhattacharjee 
for fruitful discussions and Doug Robinson for technical 
assistant during the PDF measurements. Research at the 
U. of Toronto was supported by the NSERC, CFI, and 
OMRI. This research benefited from the RIXS collabo- 
ration supported by the Computational Materials and 
Chemical Sciences Network (CMCSN) program of the 
Division of Materials Science and Engineering, U.S. De- 
partment of Energy, Grant No. DE- SC0007091. Use 
of the APS was supported by the U. S. DOE, Office of 
Science, Office of BES, under Contract No. W-31-109- 
ENG-38. Work performed at Brookhaven National Lab- 
oratory was supported by DOE, Office of Science, Divi- 
sion of Materials Science under contract No. DE-AC02- 
98CH10886. Y.-J.K. was supported by the KOFST 
through the Brainpool program. H-S Kim and J. Yu 
were supported by the NRF through the ARP (R17-2008- 
033-01000-0). H-S Kim would like to acknowledge the 
support from KISTI supercomputing center through the 
strategic support program for the supercomputing appli- 
cation research (No. KSC-2010-S00-0005). S. Manni ac- 
knowledges support from the Erasmus Mundus Eurindia 
Project. 



[1] 

[2] 
[3] 

K 

[5] 
[6] 

[r; 

[8 

[9 
[10 

[11 

[12 
[13 

[14 
[15 
[16 
[17 

[18 

[19 
[20 
[21 
[22 
[23 



Electronic address: yjkim@physic sTutoronto.cal 

Y. Okamoto, M. Nohara, H. Amga-Katori, and H. Tak- 

agi, Phys. Rev. Lett. 99, 137207 (2007). 

B. J. Kim et al., Phys. Rev. Lett. 101, 076402 (2008). 

B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, 

H. Takagi, and T. Arima, Science 323, 1329 (2009). 

G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102, 
017205 (2009). 

A. Shitade et al., Phys. Rev. Lett. 102, 256403 (2009). 
D. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 
J. Chaloupka, G. Jackeli, and G. Khaliullin, Phys. Rev. 
Lett. 105, 027204 (2010). 

J. Kim et al., Phys. Rev. Lett. 108, 177003 (2012). 
J. P. Clancy et al., e-print larXlv:1205.6540fr l (2012). 
Y. Singh and P. Gegenwart, Phys. Rev. B 82, 064412 
(2010). 

Y. Singh et al., Phys. Rev. Lett. 108, 127203 (2012). 
J. Reuther, R. Thomale, and S. Trebst, Phys. Rev. B 84, 
100406(R) (2011). 

I. I. Mazin, O. Jeschkle, R. Valneti, and D. I. Khomskii, 
e-print laiXiv: 1205.0434^ 2 (2012). 

C. H. Kim et al., Phys. Rev. Lett. 108, 106401 (2012). 
S. K. Choi et al., Phys. Rev. Lett. 108, 127204 (2012). 

F. Ye et al., Phys. Rev. B 85, 180403 (2012). 

S. Bhattacharjee, S.-S. Lee, and Y. B. Kim, New J. Phys. 
14, 073015 (2012). 

S. W. Lovesey and A. N. Dobrynin, J. Phys.: Condens. 
Matter 24, 382201 (2012). 

V. M. Katukuri et al, P hys. Rev. B 85, 220402 (2012). 
R. Comin et al., e- print larltiv:1204.447lV 2 (2012). 
M. M. Sala et al, New J. Phys. 13, 043026 (2011). 

G. Ghiringhelli et al., Phys. Rev. B 73, 035111 (2006). 
A. Uldry, F. Vernay, and B. Delley, Phys. Rev. B 85, 



125133 (2012). 

[24] G. Ghiringhelli et al., Phys. Rev. Lett. 92, 117406 (2004). 
[25] L. J. P. Ament et al., Rev. Mod. Phys. 83, 705 (2011). 
[26] L. Braicovich et al., Phys. Rev. Lett. 102, 167401 (2009). 
[27] H. Gretarsson et al., Phys. Rev. B 84, 125135 (2011). 
[28] T. Helgaker, P. j0rgensen, and J. Olsen, Molecular 

Electronic-Structure Theory (Wiley, Chichester, 2000). 
[29] M. J. OMalley, H. Verweij, and P. M. Woodward, J. Solid 

State Chem. 181, 1803 (2008), ISSN 0022-4596. 
[30] H. Kobayashi et al., J. Mater. Chem. 13, 957 (2003). 
[31] X. Liu et al., Phys. Rev. B 83, 220403 (2011). 
[32] D. S. Ellis et al., Phys. Rev. B 77, 060501 (2008). 



SUPPLEMENTAL MATERIAL 
SI. First principle electronic structure calculations 

In order to understand the origin of the low en- 
ergy excitations around 0.4 eV, we have carried out 
periodic density functional theory (DFT) electronic- 
structure calculations. We have used the DFT code 
OpenMX [si[ based on the linear combination of pseudo- 
atomic orbital(LCPAO) formalism fS2{ |. the Perdew- 
Burke-Ernzerhof GGA-functional, 8x6x8 k-points grids 
within the Brillouin zone, and 300 Ry for the real- 
space grid. The SOC is treated via a fully-relativistic 
7-dependent pseudopotential in the non-collinear DFT 
formalism SJ, £5|, S3- We applied the C2/m monoclinic 
crystal structure reported in Ref. [§3] and the zigzag- 
type magnetic order suggested in Ref. [s|| for Na2lr03. 
For Li2lr03, we used the same magnetic stucture as for 
Na2lr03 and structural data from Refs. [S9l . 



Densities of states (DOS) from GGA and GGA+U [S3j 
calculations are plotted in Fig. IS 1 1 The solid line is the 
total DOS and the filled shaded area represents partial 
d-orbitals DOS. From top to bottom, the calculations 
were done with no SOC or U, with SOC, and with both 
SOC and U. As expected the 5d states are dominant 
near the Fermi level, although there exists some modest 
amount of hybridization between Ir 5d and O 2p states. 
A moderate size U can open a (Mott) gap of 300-400 
meV, in accordance with the experimental observation. 
Compared to Na2lr03, the bandwidths in Li2lr03 are a 
little bit larger but the nature of the gap is essentially 
the same. 



S2. Quantum chemical calculations 

To investigate in detail the electronic structure and 
the essential interactions in the ^LOa iridates, we 
further performed correlated quantum chemical calcu- 
lations. In the quantum-chemical study, we employed 
multiconfiguration self-consistent-ficld and multirefer- 
ence configuration-interaction (MRCI) met hods [S 1 0( ] as 
implemented in the MOLPRO package |Sllj . The calcu- 
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FIG. SI: (color online) The evolution of the DOS for both 
Na2ir03 and Li2ir03 upon inclusion of spin-orbit coupling 
and on-site correlation. The solid line is the total DOS and 
the filled shaded area represents partial d-orbital DOS. 



lations were performed on fragments consisting of one 
central IrOg octahedron for which the local d-d transi- 
tions arc explicitly computed plus all nearest-neighbor 
(NN) octahedra, three in ^IrO^, and adjacent Na or Li 
ions. To simplify the analysis of the wave functions, the 
NN I r 4+ i ons were modeled as closed-shell Pt 4+ species 
[S12l IS 13| . The remaining part of the crystal is repre- 
sented as an array of point charges that reproduce the 
Madelung field in the cluster region. Effective core poten- 
tials and basis sets as described in earlier in vestigation s 
on Sr2±r04, Ba2ir04, and CalrOs were used |S12| . |S13| . 
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FIG. S2: Momentum dependence of the low energy RIXS 
spectra in Na2lr03 along the Q = (0 L) direction. All data 
sets collected at T — 9 K. 
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S3. Additional momentum dependence 



In many layered systems, such as the cuprates [S14j . 
momentum dependence along the L-direction is expected 
to be small. However, recent DFT calculations have 
show n dispersion of the optical gap along the L-direction 
S15| . In order to investigate this we measured the mo- 
mentum dependence along the Q = (0 L) direction 
in Na2lr03 at T = 9 K. These spectra were taken with 
the same high-resolution setup as the one in Fig. 3 (a). 
The RIXS spectra in Fig. IS2lshow no observable disper- 
sion, supporting the 2D nature of Na2ir03. We have also 
measured the momentum dependence of the high energy 
excitation (labelled D and E in Fig. 1). These spectra 
were taken with an overall resolution of ~T50 meV. Fig. 
IS3I show RIXS spectra taken along the Q = (-H -H 6.9) 
direction in Na2lr03 at T = 9 K. No observable changes 
were seen. 



FIG. S3: Momentum dependence of the high energy RIXS 
spectra in Na2lr03 along the Q = (H H 6.9) direction. All 
data sets collected at T = 9 K with a resolution of ~150 meV 
(FWHM). 



S4. X-Ray Atomic Pair Distribution Function 
(PDF) Measurements 

Atomic PDF measurements were performed at 300 K 
at the 6-ID-D beamline of the Advanced Photon Source 
at Argonne National Laboratory, utilizing a General 
Electric amorphous silicon image plate (IP) detector. A 
monochromatic incident x-ray beam 0.5mm x 0.5mm in 
size was used, conditioned to have energy of 74.353 keV 
(A=0.1668 A). Finely pulverized Li 2 Ir0 3 and Na 2 Ir0 3 
samples were packed in cylindrical polyimide capillaries 
1.0 mm in diameter and sealed at both ends. The IP de- 
tector was mounted orthogonally to the beam path with a 
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FIG. S4: (color online) Atomic PDF of Li 2 Ir0 3 and Na 2 Ir0 3 
at 300 K. (a) Li2lr03 data (gray open symbols) and refined 
C2/m model (red solid line), with the difference curve (green 
solid line) that is offset for clarity, (b) The same as (a) but 
for Na2lr03. (c) Direct comparison of Li2lr03 and Na2ir03 
experimental PDF data at 300 K, scaled such that the inten- 
sities of the first PDF peaks overlap. Dashed vertical lines 
in (a) and (b) correspond to the first two sharp PDF peaks 
in Li 2 Ir03. Arrows denote specific PDF features discussed in 
the text. 



sample to detector distance of 221.87 m m, as calibrated 
using a ceria (Ce02) standard sample [SlfJ. The 2D 
diffraction data were integrated and co nvert ed to inten- 
sity versus 29 using the software FIT2D [S17l | , where 20 is 
the angle between the incident and scattered x-ray b eam . 
The intensity data were correc ted a nd normalized |S18| 
using the program PDFgetX2 [S19j to obtain the total 
scattering structure function, F(Q) and its Sine Fourier 
transform, i.e. the atomic PDF, G(r). Q max of 25 A 
was used in the transform. The PDF analysis was carried 
out using the program PDFgui |S20l | . Atomic PDF yields 
a histogram of interatomic distances in a material, and 
provides structural informa tion on short, intermediate, 
and long range lengthscales |S18|. 

The PDF data of Li2lr03 could be explained over a 
broad r-range within the C'2/m model with all crystal- 
lographic sites fully occupied, with no indication of local 
structural distortions being present, as evident in Fig. 
[S4l (a). Structural parameters, as refined over a 1.25- 
20.0 A range, can be summarized as follows: a=5. 172(1) 
A, b=8.926(2) A, c=5. 122(2) A, = 109.91(4)°, with Ir 
at 4g (0.5, 0.167(1), 0), Lil at 2a (0,0,0), Li2 at 2d (0.5, 
0, 0.5), Li3 at 4h (0.5, 0.31(6), 0.5), 01 at 8j (0.752(5), 
0.173(2), 0.769(4)), and 02 at 4i (0.707(8), 0, 0.260(7)) 
atomic positions. The Ir-0 distances and Ir-O-Ir bond 
angles range from 2.01 to 2.04 A and 93.6 to 94.8° re- 
spectively, while Ir-Ir nearest neighbor distances are in 
the range from 2.984 to 2.992 A. 

On the other hand, while the fit of the C2/m crystallo- 
graphic model to the PDF data of Na2lr03 is consistent 



with the low-r PDF features, significant discrepancies are 
observed beginning at around 6 A, as can be seen in Fig. 
S4 (b), and no convergence could be achieved in broad r- 
rangc fits irrespective of refinement strategies attempted. 
The C2/m model is inadequate in describing the inter- 
mediate structure of Na2lr03 and the actual symmetry 
of the intermediate range structure is lower. This PDF 
is consistent with there being an appreciable amount of 
disorder in the Na2lr03 sample compared to Li 2 Ir03. To 
qualitatively illustrate this, it is useful to directly com- 
pare the experimental PDFs of the two ^IrOs systems, 
Fig. [S4](c). Local structural features are compared first. 
The first sharp PDF feature containing nearest neighbor 
Ir-0 distance in Na2lr03 occurs at around 2 A, and coin- 
cides with that of Li2lr03. The next sharp PDF feature, 
which includes the Ir-Ir near neighbor peak, in Na2lr03 
occurs at an observably larger distance in Na2lr03 (3.105 
A) than in Li2ir03 (2.985 A), indicating that the Ir-rings 
of the honeycomb network are appreciably larger, and re- 
flecting the lattice expansion on going from smaller Li to 
larger Na. Despite the shift in peak positions, the rel- 
ative intensities of these short range PDF features are 
comparable for the two samples, and have comparable 
peak widths, suggesting that the underlying local envi- 
ronments are very similar in the two systems. 

We now consider the intermediate lengthscalc in PDFs 
shown in Fig. S4 (c). What is immediately apparent is 
that starting from about 5-5.5 A the PDF features of 
Na2lr03 are visibly broader than those in Li2lr03 in- 
dicative of disorder in the former. Notably, while the 
C2/m model explains the features in the PDF of Li2ir03 
in 6-7 A range (marked by black arrows in Fig. [S4] (a)), 
the same model fails to explain the corresponding fea- 
tures in the Na2lr03 data. The two features marked by 
arrows in Li2lr03 PDF are both sharp and of approx- 
imately equal intensity. In contrast, the corresponding 
PDF peaks in Na2ir03 PDF (marked by black arrows in 
Fig. [S4](b)) have very different intensities. The PDF fea- 
ture at around 6 A in Na2lr03 is rather broad, with extra 
intensity appearing at the low-r side of the peak (indi- 
cated by red arrow in Fig. [S4] (b)), and is clearly multi- 
component. This is precisely the region where the broad 
r-range fits first fail to explain the PDF of Na2lr03, as 
these features are not incorporated into the C2 jm struc- 
tural model. Although extracting the details of the ac- 
tual intermediate structure of Na2lr03 from the present 
300 K data is a difficult task, further hampered by the 
thermal broadening that lowers the PDF resolution and 
requiring elaborate analysis of the low-temperature PDF 
data for full characterization, it is tempting to speculate 
about the origin of the observed features. Considera- 
tions of the crystallographic model reveal that in the 6 A 
range there are contributions from the third Ir-Ir neigh- 
bors in the honeycomb, as well as contributions involving 
Ir-Ir pairs belonging to two successive honeycomb sheets. 
It is plausible that the observed discrepancies originate 
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from stacking faults suggested in reference [stJ, where 
appreciable rods of diffuse scattering were observed in 
the single crystal scattering experiments. However, if 
we assume that stacking faults occur randomly at the 
< 10% level, a sharp signal in PDF such as that seen in 
the difference curve shown in Fig. [S4] (b) is difficult to 
explain. Considering that the PDF feature around 5.3 A 
in Na2lrC>3, containing the next near neighbor Ir-Ir dis- 
tance of the honeycomb, is already broad compared to its 
Li2lr03 counterpart, we cannot rule out the possibility 
that in addition to the distortions generated by stacking 
faults there arc also in-plane distortions of the honey- 
comb lattice of Na2lr03 involving correlations beyond 
that of the Ir-Ir near neighbor. 
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